Emergence of Fermi-Dirac Thermalization in the Quantum Computer Core 
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We model an isolated quantum computer as a two-dimensional lattice of qubits (spin halves) with 
fluctuations in individual qubit energies and residual short-range inter-qubit couplings. In the limit 
when fluctuations and couplings are small compared to the one-qubit energy spacing, the spectrum 
has a band structure and we study the quantum computer core (central band) with the highest 
density of states. Above a critical inter-qubit coupling strength, quantum chaos sets in, leading to 
quantum ergodicity of eigenstates in an isolated quantum computer. The onset of chaos results in 
the interaction induced dynamical thermalization and the occupation numbers well described by the 
Fermi-Dirac distribution. This thermalization destroys the noninteracting qubit structure and sets 
serious requirements for the quantum computer operability. 

PACS numbers: 03.67.Lx, 05.45.Mt, 24.10.Cn 



I. INTRODUCTION 

The key ingredient of a quantum computer WM is that 
it can simultaneously follow all of the computation paths 
corresponding to the distinct classical inputs and produce 
a final state which depends on the interference of these 
paths. As a result, some computational tasks can be 
performed much more efficiently than on a classical com- 
puter. Shor Pi constructed a quantum algorithm which 
performs large number factorization into prime factors 
exponentially faster than any known classical algorithm. 
It was also shown by Grover Q that the search of an 
item in an unstructured list can be done with a square 
root speedup over any classical algorithm. These results 
motivated a great body of experimental proposals for a 
construction of a realistic quantum computer (see ElO] 
and references therein). At present, quantum gates were 
realized with cold ions ||] and the Grover algorithm was 
performed for three qubits made from nuclear spins in 
a molecule B. For a proper operability, it is essential 
for the quantum computer to remain coherent during 
the computation process. Hence, a serious obstacle to 
its physical realization is the quantum decoherence due 
to the coupling with the external world [0-p[. In spite 
of that, in certain physical proposals, for example nu- 
clear spins in two-dimensional semiconductor structures, 
the decoherence time can be many orders of magnitude 
larger than the time required for the gate operations (see 
for example Refs. |10| , pd| ). As a result, one can analyze 
the operation of an isolated quantum computer decou- 
pled from the external world. 

However, even if the quantum computer is isolated 
from the external world and the decoherence time is infi- 
nite, a proper operability of the computer is not guaran- 
teed. As a matter of fact, one has to face a many-body 



problem for a system of n interacting qubits (two level 
systems): any computer operation -a unitary transfor- 
mation in the Hilbert space of size Nh = 2"- can be de- 
composed into two-qubit gates such as controUed-NOT 
and single qubit rotations |0,0]. Due to the unavoidable 
presence of imperfections, the spacing between the two 
states of each qubit fluctuates in some detuning interval 
6. Also, some residual interaction J between qubits nec- 
essarily remains when the two-qubit coupling is used to 
operate the gates. 

In [|l2|-|l4| an isolated quantum computer was modeled 
as a qubit lattice with fluctuations in individual qubit 
energies and residual short-range inter-qubit couplings. 
Similarly to previous studies of interacting many-body 
systems such as nuclei, complex atoms, quantum dots, 
and quantum spin glasses |15|-pi[ , the interaction leads to 
quantum chaos characterized by ergodicity of the eigen- 
states and level spacing statistics as in Random Matrix 
Theory 123] . The transition to chaos takes place when the 
interaction strength is of the order of the energy spacing 
between directly coupled states ||l^,|l^-^. This border 
is exponentially larger than the energy level spacing in a 
quantum computer. 

In this paper we show that the onset of chaos leads to 
occupation number statistics given by the Fermi-Dirac 
distribution. This means that a strong enough interac- 
tion plays the role of a heat bath, thus leading to dy- 
namical thermalization for an isolated system. In such 
a regime, a quantum computer eigenstate is composed 
by an exponentially large (with n) number of noninter- 
acting multi-qubit states representing the quantum regis- 
ter states. As a result, exponentially many states of the 
computation basis are mixed after a chaotic time scale 
and the computer operability is destroyed. We note that 
the dynamical thermalization has been discussed in other 



many-body interacting systems in [ [16| , |17[|1S|J23| . 

The paper is composed as follows. In Section |i we de- 
scribe our qubit lattice quantum computer model [ 12|-p^ ; 



in Section [II we discuss the statistical properties of the 
eigenvalues of this model; in Section |^ we study the oc- 
cupation number distribution and compare different def- 
initions for the effective temperature of the system; the 
conclusions are presented in Section M. 



II. THE MODEL 



We consider a model of n qubits on a two-dimensional 
lattice with nearest neighbors intcr-qubit coupling P4| . 
The Hamiltonian of this model, introduced in ||l^, reads: 



H 



Z^ 



T,af 






(1) 



where the ai are the Pauli matrices for the qubit i and 
the second sum runs over nearest-neighbor qubit pairs on 
a two-dimensional lattice with periodic boundary con- 
ditions applied. The energy spacing between the two 
states of a qubit is determined by Ti = Aq -I- 5i, with 
6i randomly and uniformly distributed in the interval 
[—6/2,6/2]. Therefore the detuning parameter 6 gives 
the width of the Ti distribution around its average value 
Aq. For generality we choose the couplings J^, which 
represent the residual interaction, randomly and uni- 
formly distributed in the interval [—J, J]. The model 
(|l|) can be considered as a standard generic quantum 
computer model, in which the unavoidable system im- 
perfections generate a residual inter-qubit coupling and 
energy fluctuations. In a sense (|l|) describes the quantum 
computer hardware, while to study the gate operations in 
time one should include additional time-dependent terms 
in the Hamiltonian. At J = 0, the noninteracting eigen- 
states of the model can be written as {ipk) = |q^i, •■■, Q^n), 
where ai — 0,1 marks the polarization of each individ- 
ual qubit. These are the ideal multi-qubit eigenstates of 
a quantum computer, the quantum register states used 
for computer operations. For J ^ 0, these states are no 
longer eigenstates of the Hamiltonian, and the new multi- 
qubit eigenstates are now linear combinations of different 
quantum register states. 

Here we focus on the case 6 ^ Aq, which corresponds 
to the situation where fluctuations induced by imperfec- 
tions are relatively weak p^ . In this case, the unper- 
turbed energy spectrum of (|lj) (corresponding to J = 0) 
is composed of n -I- 1 well separated bands, with inter- 
band spacing 2Ao. Since the 6i randomly fluctuate in an 
interval of size 6, each band at J = (except the extreme 
ones) has a Gaussian shape of width « y^6 [E5[. The 
average number of states inside a band Nb is of the order 
of Nh /n = 2"/n, so that the energy spacing between ad- 
jacent multi-qubit states inside one band is exponentially 
smaU: A„ ~ n^/^2-'^6. 



In the presence of a residual interaction J ^ 6, the 
spectrum still has the above band structure with an ex- 
ponentially large density of states. For J, (5 <C Aq, the 
interband coupling is very weak and can be neglected. In 
this situation, the Hamiltonian (n^) is, to a good approx- 
imation, described by the renormalized Hamiltonian 



n+l 



Hp = Y.PkHPk, 



(2) 



fc=i 



where Pk is the projector on the fc*'* band, so that qubits 
are coupled only inside one band. We concentrate our 
studies on the central band. For an even n this band 
is centered exactly at E = 0, while for odd n there are 
two bands centered at -E = ±Ao, and we consider the 
one a.t E = — Aq. The central band corresponds to the 
highest density of states, and in a sense represents the 
quantum computer core: an exponentially large number 
of states allows to take advantage of quantum parallelism 
in computer algorithms ||U-0]. On the other hand, quan- 
tum chaos and ergodicity first appear in this band, which 
therefore sets the limit for operability of the quantum 
computer. Inside this band, the system is described by 
the renormalized Hamiltonian Hp which depends only on 
the number of qubits n and the dimensionless coupling 
J/6. 



III. SPECTRAL STATISTICS 

The results of Refs. pslfll showed that the quantum 
chaos border in (nl) corresponds to a critical interaction 
Jc given by: 



Jc 



C6 
n 



(3) 



where C is some numerical constant. This border is expo- 
nentially larger than the energy spacing between multi- 
qubit states A„ . This is in agreement with previous stud- 
ies of complex interacting many-body systems |l5| , [18| -pi[ , 
in which the transition to quantum chaos takes place 
when the interaction matrix elements between directly 
coupled states become larger than their energy spacing. 
Since the interaction is of a two-body nature, each nonin- 
teracting multi-qubit state \ipk) has nonzero coupling ma- 
trix elements only with about n other multi-qubit states. 
Therefore, the number of directly coupled states is much 
smaller than the number of multi-qubit states inside the 
central band, Np — n\/{[n/2]\{n — [n/2])!) (we consider 
the band with the number of spins up given by the inte- 
ger part of n/2). These couplings induce transitions in 
an energy interval of order 5 (we assume that J is of the 
order of or smaller than 6). Therefore the energy spacing 
between directly coupled states is Ac ~ 6 /n. The tran- 
sition to chaos takes place for J — Jc^^ /S.c, which leads 
to the relation (^. 
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FIG. 1. Dependence of r; on the scaled coupling Jn/5, 
for n — 9 qubits (circles, Nd = 10* random realizations of 
Si,Jij), n = 12 (squares, Nn ~ W^), n — 15 (diamonds, 
Nd = 45), n — 16 (empty triangles, Nd = 23), and n = 18 
(filled triangles, Nd = 3). The parameter rj is computed for 
±5% of states around the center of the energy band. 



The transition to quantum chaos and ergodic eigen- 
states can be detected in the change of the spectral statis- 
tics of the system. A convenient way is to look at the level 
spacing statistics P{s), which gives the probability to find 
two adjacent levels whose spacing, normalized to the av- 
erage level spacing, is in [s,s + ds]. In fact, P{s) goes 
from the Poisson distribution Pp{s) = exp(— s) for non- 
ergodic states to the Wigner-Dyson distribution P\y (s) = 
(7rs/2) exp(— 7rs^/4), corresponding to Random Matrix 
Theory, for ergodic states p3]. To analyze the change of 
P{s) with the coupling J one can conveniently use the pa- 
rameter 77 = J^°{P{s)~Pw{s))ds/ J^"{Pp{s)-Pw{s))ds 
|l9| , where sq = 0.4729... is the first intersection point 
of Pp{s) and Pw{,s). In this way Pp{s) corresponds 
to 77 = 1 and Pwis) to rj — 0. Fig|| gives the depen- 
dence of the parameter 77 on the scaled coupling Jn/S 
at different system sizes, for states near the middle of 
the energy spectrum (±5% of levels around the band 
center). In order to reduce statistical fluctuations, we 
use 3 < Np) < 10^ random realizations of Si,Jij. In 
this way the total number of spacings Ns is varied in 
the interval 1.4 x 10^ < iV^ < 1.2 x 10^ The Pois- 
son to Wigner-Dyson crossover becomes sharper when 
n increases, suggesting a sharp transition in the ther- 
modynamic limit. We note that, since the number of 
random 6i, Jij values is not large, significant fluctuations 
are present in the 77 curves when one changes the ran- 
dom realization. For example, for 77 = 18, at J = 0.26, 
one has 77 = 0.44, 0.29, 0.06 in the Np) = 3 random re- 
alizations considered. This reflects the general property 
of fluctuations which become stronger near the critical 
transition point. On the other hand, even if the number 
of considered spacings is always large, for 77 = 15, 16 and 



especially for n = 18, we have only a small number of 
random realizations, due to the very slow convergence of 
the Lanczos algorithm p6[ near the band center, where 
the density of states becomes exponentially large. These 
fluctuations prevent us from precisely evaluating the crit- 
ical scaled coupling for the Poisson to Wigner transition. 
The minimum spreading of curves is for 77 (Jc) « 0.2, cor- 
responding to Jcn/6 ~ 3.7, in agreement with the results 
of Rcf . Q . We stress that the chaos border is exponen- 
tially larger than the multi-qubit level spacing, e.g., for 
n = 18, Jc ~ 0.26 > A„ « 7 X 10-^(5. 
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FIG. 2. Level spacing statistics corresponding to Fig. hi at 
n = 16, J = 0.055 (circles, rj = 1.01), J = 0.25 (triangles, 
rj = 0.32), and J = 0.45 (squares, r/ = 0.05). Full curves show 
Poisson and Wigner-Dyson distributions. 



The level spacing statistics near the band center is 
shown in FigJ| at different coupling strengths J for 
n = 16. The transition from the Poisson to the Wigner- 
Dyson statistics is evident. In Fig.|| we show that, even 
away from the band center, the transition happens at ap- 
proximately the same critical coupling. This is due to the 
fact that, even if the 77-body density grows exponentially 
with the excitation energy, the density of coupled states 
remains roughly the same, except near the band edges 
pTf . We note that the data are shown only for half of 
the central band {E < 0), since the density of states is 
symmetric around E = due to the presence of an upper 
bound in the single qubit energies. 
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FIG. 3. Dependence of rj on the energy E (normalized to 
the total band width B taken from numerical data) for n — 16, 
J = 0.055 (circles), J = 0.155 (squares), J — 0.2S (diamonds), 
and J — 0.45 (triangles). Data are averaged over Nd ~ 15 
random realizations. 



IV. DYNAMICAL THERMALIZATION 

The transition in the level spacing statistics reflects 
a qualitative change in the structure of the eigenstates 
[^,|3[. While for J ^ Jc the eigenstates are very close 
to the quantum register states, for J > Jc each eigenstate 
1 0m) becomes a superposition of an exponentially large 
number of noninteracting eigenstates \ipk)- It is conve- 
nient to characterize the complexity of an eigenstate \4>rn) 
by the quantum eigenstate entropy 



Nb 

^ Wk,n l0g2 Wkr. 

k=l 



(4) 



where Wkm is the probability to find the quantum reg- 
ister state \ipk) in the eigenstate \(j)m) of the Hamilto- 
nian (Wkm = |(^fc|(/>m)P)- In this way 5g = if \(j)„i) 
is a quantum register state (J = 0), 5'g = 1 if 10m) is 
equally composed of two \ipi), and the maximum value is 
Sq — log2 Nb if all states equally contribute to |0m)- 

Above the chaos border (J > Jc) one eigenstate is 
composed of order 2^1 quantum register states, mixed 
inside the Breit-Wigner width T ^ J^/Ac ~ J^n/S [^. 
As a result, the residual interaction disintegrates a quan- 
tum register state over an exponentially large number of 
states after a chaotic time scale r^ ~ l/F |0,p8[. After 
this time the quantum computer operability is certainly 
destroyed, unless one can apply quantum error-correcting 
codes (see Mpl and references therein) operating on a 
shorter time scale. This destruction takes place in an 
isolated system, without any external decoherence pro- 
cess. It happens due to inter-qubit coupling, which can 
mimic the effect of a coupling with the external world. In 
the following we show that in the quantum chaos regime 



a statistical description of our isolated n-qubit system is 
indeed possible, similarly to results found for other phys- 
ical systems in ||g 17,2^]. 

We concentrate on the distribution of the occupation 
numbers rii, defined as the probability that the qubit 
(spin) at the site i is in its up polarization state. Given 
an eigenfunction |0m) with eigenvalue Em, one can write: 



Nb 



i(m) = ^M^fcm(V'fcl^ilV'fc), 



(5) 



k=l 



where hi is the occupation number operator, and the 
term (V'fel^ilV'fe) equals 1 or depending on whether the 
spin at the site i is up or down. 

For noninteracting qubits one can write, e.g. for the 
central band. 



1=1 

n 

^n,(m)(5,; =-B„, (nj(m)=0, 1), 



(6) 



where < = Sm/2 + Y.^ <5^/2 (^m = E.(2"^(fc) - l)5^). 
As ni{ra) = 0, 1, the relations (^ are the usual ones used 
to derive the Fermi-Dirac distribution for an ideal gas of 
many noninteracting particles in contact with a thermo- 
stat. However, here we consider an isolated system of 
relatively few interacting particles. Nevertheless, recent 
studies ||l^,|l^,^ have demonstrated that interaction can 
play the role of a heat bath, thus allowing one to use a 
statistical description even in an isolated system with few 
particles. The Fermi-Dirac statistics appears due to the 
fact that the number of spins up/down is fixed and in 
this way they become equivalent, for the purposes of a 
statistical description, to electrons/holes. 

In Fig.y we show the occupation number distribution, 
averaged over a few consecutive levels |0m) and over 
iVxj = 100 random realization of 5^, J^j, for n = 16 qubits, 
both in the integrable regime (top figures) and in the 
quantum chaos regime (bottom figures). We see that 
in both cases this averaged distribution is in very good 
agreement with the Fermi-Dirac distribution 



FD 



exp(/?((5, + <5/2-m)) + 1' 



(7) 



where /i is the chemical potential and j3 = l/Tpu is the 
inverse temperature (we set fcs — 1). Taking into ac- 
count the constraint set by the fixed number of spins up 
(Yl^nf^ = [n/2]), TpD is the only fitting parameter. 

The goodness of the fit given by the expression (R) is 
hardly surprising as statistical distributions are obtained 
for noninteracting particles after a correct counting of 
states. In this procedure, a weak interaction gives a slight 
increase of the system temperature p7| . 



A good agreement between the numerical data in Fig.^ 
and the theoretical distribution (Q) does not mean that 
automatically there is equilibrium and thermalization for 
a given realization. This is outlined in Fig.@, which 
shows the occupation numbers for a single eigenstate 
of a given random realization. In the upper figures 
(J = 0.03(5 << Jc ~ 0.2(5) a given eigenstate signifi- 
cantly projects only over a single quantum register state 
{Sq < 1) and therefore half of the occupation numbers 
is close to 1, half close to 0, and the Fermi-Dirac distri- 
bution (M) is very far from the actual distribution. On 
the contrary, in the quantum chaos regime (lower figures, 
J = 0.3(5 > Jc), where a large number of quantum regis- 
ter states are mixed in a single eigenstate, there is a good 
agreement between the occupation number distribution 
and the Fermi-Dirac distribution. 
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FIG. 4. Distribution of the occupation numbers rii as a 
function of the qubit detunings Si, for n — 16 qubits, cou- 
pling strength J = 0.03(5 (upper figures) and J = 0.3(5 
(lower figures), levels m = 5 — 10 (left) and m = 95 — 100 
(right). Solid lines give Fermi-Dirac fits with effective temper- 
ature Tfd- Top left: Tfd = 0.13(5, mean excitation energy 
SE = 0.215, mean quantum eigenstate entropy Sq = 0.56; 
top right: Tfd = 0.20(5, SE = 1.145, 5", = 0.81; bottom 



left: Tfd = 0.175, SE = 0.245, S„ 



6.03; bottom right: 



Tfd = 0.255, 5E = 1.325, 5, = 8.04. Data are averaged over 
Nd = 100 random realizations. 

In order to make quantitative the comparison with 
the Fermi-Dirac distribution, we introduce a parameter 
which measures the root mean square deviation of the 
actual distribution from (^: 



(7FD{m) 



\ 



1 " 

— 2j('^i("i) ~ nf^{m)Y. 



(8) 



For the case of Fig.H, we have apo — 6.3 x 10~^ (top 
left), apD = 8.8 X 10"^ (top right) and much lower values 
above the chaos border: apD = 1-4 x 10^^ (bottom left), 



o'FD — 1-7 X 10~^ (bottom right). The maximum value 
fj-rna^ _ Q_g jg obtained at the band center {Tpu = oo) 
for J — 0, when Ui = 1 for [n/2\ spins and rii — for the 
remaining ones. 

We introduce the thermalization border Jt as follows: 
for J < Jt eigenstates close in energy yield completely 
different Tii-distributions, for J > Jt the n.^-distribution 
is stable with respect to the choice of a specific eigen- 
state in a small energy window. The appropriate quan- 
tity to be considered in addition to apD is the root mean 
square deviation of the occupation numbers for consecu- 
tive eigenstates: 



\ 



1 " 

- y^{ni{m + 1) ~ ni{m)y 

n ^ — ^ 



(9) 



For the case of Fig.^, we have Ug — 0.1 (top left), 
as — 0.15 (top right) and much lower values in the chaotic 
regime: as = 2.9 x lO'^ (bottom left), as = 3.3 x 10"^ 
(bottom right). 

The conclusions drawn from Fig.o are also applied to 
Fig.p where, thanks to the effectiveness of the Lanczos 



algorithm |26| near the band edges, it was possible to 
consider a larger number of spins (n = 24, corresponding 
to a very large Hilbert space dimension Nh ~ 1.7 x 10^, 
with Nb ~ 2.7 X 10^ levels in the central band). 
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FIG. 5. Same as in Fig.Wbut for a given random realization 
and a single eigenstate for n = 16 qubits (left: m — 5; right: 
m = 100). Top left: J = 0.035, Tfd = 0.085, SE = 0.255, 
Sg = 0.23; top right: J = 0.035, Tfd = 0.155, SE = 0.975, 

0.35, Tfd = 0.095, SE = 0.285, 
0.205, SE^ 1.195, 



Sq = 0.49; bottom left: J 

Sg = 5.85; bottom right: J = 0.35, Tfd 

Sa = 8.41. 



Fig.M shows the parameters apo and as as a func- 
tion of the energy, for different values of the coupling 
strength J. As for the transition to chaos (Fig.m, the 
thermalization occurs at approximately the same criti- 
cal coupling also away from the band center. We note 



that our data give CTs ~ V^f^FD, which can be eas- 
ily understood as follows: X]i("i(™ + 1) ~ 'ni{m))^ « 
E,(n,(m)-nf^(m))2 + ^^(n,(m + l)-nf^(m+l))2- 
2 J2iM'm')-n[^{m)){n,{m+l) - nf ^(m+ 1)), the last 
term in the sum becoming negligible after ensemble aver- 
aging. We stress that the very good agreement between 
the parameters apD and CTs implies that the Fermi-Dirac 
distribution emerges when the occupation number dis- 
tribution is statistically stable, i.e. the system is ther- 
malized. The curve for J = 0.05(5 lowers near the band 
edge since for a small excitation energy 6E a small num- 
ber rieff ~ {nSE/Sy^'^ of spins ("fermions") is available 
for fluctuations in the vicinity of the Fermi level (see the 
note 13). 
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FIG. 6. Same as in Fig.g but for n = 2A qubits. Top left 
J = 0.055, Tfd = 0.055, SE = 0.16(5, Sq = 0.51; top right 
J = 0.055, Tfd = 0.125, SE = 0.645, 5"^ = 1.84; bottom left 



J 
J 



- 0.45, Tfd 
■ 0.45, Tfd 



0.065, SE : 
0.195, SE 



0.385, 5", 
1.165, So 



7.75; bottom right 
12.55. 



It is interesting to compare the temperature Tpu ob- 
tained from the Fermi-Dirac fit with different definitions 
of temperature ||l^,|l^,^ , which are known to be equiv- 
alent at the thermodynamic limit. First of all we use the 
canonical expression 



E{Tcan) — 



_i \ ca7i / 



Ne 



^exp 



J- cat 



(10) 



where Em are the exact eigenenergies of the interact- 
ing system. The very good agreement between Tp^ and 
Tcan (see Fig.ph supports the validity of a statistical de- 
scription for our isolated quantum computer model. This 
means that in such closed system the inter-qubit resid- 
ual interaction plays the role of a heat bath in an open 
system. In particular, we expect that residual interac- 



tion mimics to a certain extent the effect of coupling to 
external world. 
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FIG. 7. Root mean square deviation ctfd of the occupation 
number distribution with respect to the Fermi-Dirac fit as a 
function of the energy E (normalized to the band width B), 
for n = 16, J — 0.055 (circles), J = 0.155 (squares), J — 0.25 
(diamonds), and J — 0.45 (triangles). Lines give (Ts/\/2, with 
(Js root mean square deviation of the occupation numbers for 
consecutive levels. Data are averaged over Nd ~ 2 random 
realizations. 
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FIG. 8. Dependence of different definitions of temperature 
T on the scaled energy E/B, for n = 16, J = 0.35, Nd = 2: 
Tfd (circles), Tcan (full curve), and Tth (dashed curve). 

Finally we compare the effective temperature Tpu with 
the standard thermodynamic temperature Tth, defined 
by 



1 

Tth 



dS, 



th 



dE' 



dlnp 
dE' ' 



(11) 



where Sth = Inp is the thermodynamic entropy, with p 
density of states (here £" replaces S^). A Gaussian fit, 



p = (l/-\/27rcr)exp(— £"^/cr^), provides an excellent ap- 
proximation to the actual data for the density of states, 
with the exception of the first few levels near the band 
edges. Therefore one gets Tth = — cr^/£", with cr^ vari- 
ance of the Gaussian fit. Fig.@ shows that the difference 
between Tth and Tpo ~ Tcan increases as one moves away 
from the band center. This is due to the fact that, con- 
trary to the noninteracting case (see Eq.(^)), the actual 
system energy E' is different from ^^ riiSi, since the total 
energy is renormalized due to interaction p7[ . We remark 
that, not only on average but also for a single eigenstate 
of a given realization, the agreement between Tpjj and 
Tcan is good, while Ttn becomes closer to Tpu and Tcan 
when the excitation energy increases. In the case of Fig.^ 
we obtain Tpo = 0.08(5, Tcan = 0.07(5, Tth = 0.19(5 (top 
left), Tfd = 0.15(5, Tcan = 0.16(5, Tth = 0.24,5 (top right), 
Tfd = 0.09(5, Tcan = 0.08(5, Tth = 0.22<5 (bottom left), 
Tfd = 0.20(5, Tcan = 0.19<5, Tth = 0.29(5 (bottom right). 
We also note that the effective temperature diverges at 
the band center and is negative in the upper part of the 
spectrum. This is typical of models with an upper bound 
in the single particle energies, in which the density of 
states has a maximum. 




FIG. 9. Dependence of apD on the scaled coupling Jn/5, 
for n — 9 (circles, Nd = 2x 10^ random realizations oiSi, Jij), 
n = 12 (squares, TVu = 2x 10^), n = 15 (diamonds, Nd = 20), 
and n = 16 (triangles, Nd = 10). The parameter apD is com- 
puted for ±5% of states around the band center. Transition 
to thermalization takes place at Jtn/5 « 3.2. 



The dependence of the thermalization parameter apD 
on the scaled coupling Jn/S at different system sizes 
is shown in Fig.0 (very similar results are obtained for 
(Ts/\/2). Our data show that the crossover to a thcrmal- 
ized distribution sharpens when the number of qubits 
increases, in a way consistent with a sharp thermaliza- 
tion border in the thermodynamic limit, at Jtn/5 ~ 3.2. 
The similarity between the results of Fig.H^ {Jcu/S w 3.7) 
and Fig.|| leads us to the conclusion that the chaos bor- 



der coincides with the thermalization border. This looks 
quite natural since the Poisson level spacing statistics 
indicates the existence of uncoupled parts in the whole 
system, thus preventing thermalization. On the contrary, 
in the chaotic regime each eigenfunction spreads over an 
exponentially large number N « 2^1 of quantum reg- 
ister states, resulting in the Wigner-Dyson statistics. In 
this regime the fiuctuations of eigenstate components are 
Gaussian [n^ and therefore, according to the central limit 
theorem, the fiuctuations of the occupation numbers are 
small: Arii oc N~^/^ ^ 1. For this reason eigenstates 
close in energy give similar n^-distributions, which means 
that there is equilibrium in the statistical sense. 



V. CONCLUSIONS 

The results presented in this paper show that the resid- 
ual inter-qubit coupling can lead to quantum chaos and a 
statistical description of the occupation numbers in close 
agreement with the Fermi-Dirac distribution. We stress 
that the transition to quantum chaos is an internal pro- 
cess which happens in a perfectly isolated system with no 
coupling to the external world. Nevertheless, the ther- 
malization which appears in this closed system due to 
inter-qubit coupling can mimic the effect of an exter- 
nal thermal bath. Above the chaos/thermalization bor- 
der the quantum register states are destroyed after the 
chaotic time scale and therefore serious restrictions are 
set to the quantum computer operability. However, be- 
low this border the imperfection effects only slightly dis- 
turb the ideal multi-qubit states and the quantum com- 
puter can operate properly in this regime. 

We thank Bertrand Georgeot for stimulating discus- 
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